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ABSTRACT 



We present the results of hydrodynamical simulations of the disk photosphere 
irradiated by strong X-rays produced in the inner- most part of the disk of an 
accreting black hole. As expected, the irradiation heats the photosphere and 
drives a thermal wind. To apply our results to the well-studied X-ray transient 
source GRO J1655-40, we adopted the observed mass of its black hole, and 
the observed properties of its X-ray radiation. To compare the results with the 
observations, we also computed transmitted X-ray spectra based on the wind 
solution. Our main finding is: the density of the fast-moving part of the wind is 
more than one order of magnitude lower than that inferred from the observations. 
Consequently, the model fails to predict spectra with line absorption as strong 
and as blueshifted as those observed. However, despite the thermal wind being 
weak and Compton thin, the ratio between the mass-loss rate and the mass 
accretion rate is about seven. This high ratio is insensitive to the accretion 
luminosity, in the limit of lower luminosities. Most of the mass is lost from the 
disk between 0.07 and 0.2 of the Compton radius. We discovered that beyond 
this range the wind solution is self-similar. In particular, soon after it leaves 
the disk, the wind flows at a constant angle with respect to the disk. Overall, 
the thermal winds generated in our comprehensive simulations do not match the 
wind spectra observed in GRO J1655-40. This supports the conclusion of Miller 
et al. and Kallman et al. that the wind in GRO J1655-40, and possibly other 
X-ray transients, may be driven by magnetic processes. This in turn implies that 
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the disk wind carries even more material than our simulations predict and as 
such has a very significant impact on the accretion disk structure and dynamics. 

Subject headings: accretion, accretion disks - methods: numerical - HD 

1. INTRODUCTION 

Most X-ray sources are powered by disk accretion onto compact objects. Therefore 
a main challenge for X-ray astronomy is to understand the mechanisms that enable this 
process. High quality and high spectral-resolution observations obtained with Chandra, 
XMM-Newton, and Suzaku allow us to study disk accretion and related outflows better 
than ever before. These observations are especially revealing if taken of relatively bright, 
well studied objects such as X-ray transient sources, with GRO J1655-40 being a very fine 
example. 

Many properties of GRO J1655-40 are well constrained. For example, GRO J1655-40 
is a black hole binary at a distance of 3.2 kpc containing a black hole with a mass of 7.0 
M0 that accretes from an F3 IV-F6 IV star with a mass of 2.3 Mq in a 2.6-day orbit. The 
inner disk is viewed at an inchnation of 67° - 85° (nearly edge-on; Orosz & Bailyn 1997). 
When GRO J1655-40 is in an X-ray-bright phase, it is possible to obtain with Chandra high 
signal-to- noise spectra as shown by Miller et al. (2006a; M06 hereafter). In fact, the quahty 
of the spectra obtained by M06 is good enough to reveal many absorption fines significant 
at the 5 a level of confidence or higher. Over 70 of these lines can be identified as resonance 
lines from over 32 charge states. This is in sharp contrast to high resolution spectra of other 
systems (see below) and of GRO J1655-40 in other spectral states, when only the absorption 
lines of Fe XXV and XXVI are detected (e.g.. Miller et al. 2006b; Kubota et al. 2007; Miller 
et al. 2008; Nielsen & Lee 2009). 

It is hoped that such high quality spectroscopy can provide new and surprising results. 
Indeed, absorption lines discovered by M06 are blueshifted and likely produced in a disk 
wind, making GRO J1655-40 the best case for an X-ray binary with an X-ray absorbing 
disk wind. But there are also other cases. For example, Chandra and RTXE spectroscopy 
of the microquasar H 1743—322 reveals blueshifted absorption fines that are fikely formed 
in a disk wind (Miller et al. 2006b). Another example is the X-ray transient 4U 1630-472. 
Outburst spectra of this source obtained with Susaku, show iron absorption lines indicative 
of a disk wind (Kubota et al. 2007). More recently, Neilsen & Lee (2009) found an X-ray 
absorbing disk wind in the microquasar GRS 1915+105. In addition, there is an outflowing 
X-ray absorber in Circinus X-1, which could either be a disk wind or a wind from a massive 
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companion star (Schulz & Brandt 2002). 

The disk wind discovered in GRO J1655-40 gives us a good testbed for constraining wind 
properties using X-ray observations. For example, by fitting to observations 1-dimensional 
(1-D) models, M06 and Miller et al. (2008) put strict limits on the ionization balance in 
the wind. This was facihtated by the fact that some of the detected fines provide a density 

diagnostic. The main results from the fitting yielded the following parameters: number 
density n between 5 x 10^^ cm^"^ and 2 x 10^"^ cm^'\ column density A^h = 7.4 x lO^"*^ cm~^, 
distance from the central object to the point of emerging of the strong wind 4.8 x 10^ cm 
and velocity = 500 kms~^. The lines show blueshifts in the 300-1600 km s~^ range. 

Using these wind properties, together with the observed wind speed and the system 
luminosity, M06 concluded that the inferred wind location is well within the Compton radius 
(the radius where the gravitational and thermal pressures are equal), defined by: 

GMBKmpH 

^ mc ' 

where Mbh is the mass of the black hole, /i is the mean molecular weight, rrip is the proton 
mass, and Tic is the Compton temperature (for Mbh = 7 M©, /i = 0.6, and Tic = 1-4 x 10'' K, 
the Compton radius is i?ic = 4.8 x 10^^ cm). This wind location is then inconsistent with 
an outfiow being driven by thermal expansion, even though weak fiows are possible at ~ 0.1 
of the Compton radius, as found by Woods et al. (1996; see also Proga & Kallman 2002, 
hereafter PK02). On the other hand, the wind cannot be driven by radiation pressure in 
this system because the luminosity, L^, relative to the Eddington limit, L/L^dd = 0.03, is 
too low (e.g., PK02). Therefore, by a process of elimination, M06 concluded that the only 
plausible mechanism that could drive the wind is magnetic processes. M06 also concluded 
that disk accretion itself is driven by magnetic fields. 

Although magnetic forces can drive a disk wind (e.g., Blandford & Payne 1982) and 
magnetic fields are a very important ingredient of accretion disks (Balbus & Hawley 1998), 
M06's arguments for magnetic driving are indirect (see also Proga 2006). Ideally, one should 
demonstrate that magnetic forces can drive a disk wind capable of reproducing the observed 
spectra. However, that is a relatively challenging task, and instead one could try to verify first 
whether thermal driving is indeed unsatisfactory. This is a very relevant question especially 
in light of the work done by Netzer (2006; N06 hereafter), who argued that it is possible to 
produce a simple thermally driven wind that is consistent with the observed properties of 
GRO J1655-40. 

One can argue for or against thermal driving by comparing the escape velocity and the 
isothermal sound speed at the Compton temperature (equivalently, one can compare the 
wind launching radius, Ri with, i?ic as M06 did). To use this basic physical argument, one 
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needs to estimate the wind temperature and Ri. The former can be constrained relatively 
well for the Compton heated gas, whereas the latter is much more difficult to constrain. In 
fact, the controversy about the role of thermal driving has its origin in two different estimates 
for Ri. M06 estimated Ri = 5x 10^ cm, whereas N06 estimated Ri = 5x 10^° — 5 x 10^^ cm. 
These two estimates differ by 2-3 orders of magnitude! For the lower estimate, the wind 
cannot be thermally driven, whereas for the higher estimate it can. We note that the two 
groups also found different values for Ric (i.e., M06's estimate is 7 x 10^^ cm whereas N06's 
estimate is 5 x 10^^ cm). Additionally, the two groups disagree about the gross properties 
of the thermal wind. In particular, contrary to N06, M06 claimed that both theory and 
simulations (Begelman et al. 1982; Woods et al. 1996) predict mass-loss rates of thermal 
winds that are much too low to account for the observed density and blueshift. 

The reason for the above disagreement can be traced down to the differences about 
the absorption fines of Fe XXII and Fe XXIII. N06 assumed saturated lines and the absolute 
covering factor of the X-ray source to be 0.75, whereas M06 took the lines to be unsaturated 
and used a covering factor of 1. 

Miller et al. (2008) confirmed their previous results and stated that if N06's claim of 
the Fe XXII line being saturated were correct then the ratios of some other fines should be 
different than what is observed. Miller et al. also re- iterated their point that that FeXXII 
lines arc optically thin and that the covering factor along the line of the sight is near unity. 
They also stated that N06 used too low a density (by factor of 5-10) as a consequence of his 
line saturation assumption. 

The conclusion of M06 and Miller et al. (2008) was confirmed by Kallman et al. (2009) 
who used spectral fitting to show that the ionization conditions in a 1-D model are not 
consistent with the wind being driven by thermal expansion. 

Generally, M06, NOG, and Kallman et al. (2009) agree that thermal driving operates in 
this system but disagree about the properties of the thermal wind. Overall, it appears that 
thermal driving is unfavorable. However, this conclusion is based on simplified 1-D wind 
models and needs to be confirmed by a detailed physical model that takes into account the 
intrinsic multi-dimensional geometry of disk winds. 

In this paper, we present 2.5-D axisymmctric, time-dependent hydrodynamical numer- 
ical simulations of thermally driven disk winds. We focus on the results obtained after the 
simulations reached the steady state and compare the results with observations. Our simu- 
lations do not include any magnetic processes and therefore do not address the possibility of 
the wind being magnetically driven, but they instead focus on the thermal contribution to 
the wind. The main goal of our simulations is to assess within the assumptions whether or 
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not the thermally driven wind can account for the observations. We also consider possible 
implications of the fact that the mass-loss rate of the wind is several times higher than the 
total mass- accretion rate. 

The paper is organized in the following way: Section 2 describes the methods used; 
Section 3 lays out the initial conditions of the simulation and discusses the properties of the 
fiducial run as well as its differences from other runs; Section 4 compares the results to the 
observations; Section 5 summarizes and discusses our results. 



2. METHOD 

2.1. HYDRODYNAMICS 

To compute the structure and evolution of a disk irradiated by the central source, we 
numerically solve the equations of hydrodynamics 

^ + PV.V = 0, (2) 
P^ = -VP + pg (3) 
P^J-) ^-PV-v + pC, (4) 



Dt Vp, 

where p is the mass density, P is the gas pressure, v is the velocity, e is the internal energy 
density, JC is the net cooling rate, and g is the gravitational acceleration of the central 
object. We adopt an adiabatic equation of state, P — (7 — 1) e, and consider models with 
the adiabatic index, 7 = 5/3. 

We use the ZEUS-2D code (Stone & Norman 1992) extended by Proga et al. (2000) to 
solve eqs. 2-4. We perform our calculations in spherical polar coordinates (r, 9, 0) assuming 
axial symmetry about the rotational axis of the accretion disk {9 — 0°). 

Our computational domain is defined to occupy the angular range 0° < ^ < 90° and the 
radial range ri < r < r^- The numerical resolution consists of 200 cells in the r direction 
and 100 cells in the 9 direction. In the angular direction, we used the following ratios: 
d^k+i/c?^k=0.95, 0.97, 0.99, and 1.00 [i.e., the zone spacing increases towards the pole]. 
Gridding in this manner ensures good spatial resolution close to disk. In the r direction, we 
used cirk+i/c?rk — 1.04, which enables a good resolution at smaller radii. 

We chose the boundary condition at the pole (i.e., ^ = 0°) as an axis-of-symmetry 
boundary condition. At the disk (i.e., 9 — 90°), we apphed the refiecting boundary condition. 
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For the inner and outer radial boundaries, we used an outflow boundary condition (i.e., to 
extrapolate the flow beyond the boundary, we set values of variables in the ghost zones equal 
to the values in the corresponding active zones, see Stone & Norman 1992 for more details). 



2.2. RADIATION FIELD 

We deal with the radiation field and radiation heating and cooling in the the same 
manner as described by Proga et al. (2002; see also Proga et al. 2000). 

The net coohng rate is a function of the photoionization parameter, which is defined as: 

(5) 

where J-'x is the local X-ray flux, n = p/ {m^ii) is the number density of the gas. We consider 
fully ionized gas with ji — 0.6. The local X-ray flux is corrected for the optical depth effects: 

Jk = exp(-Tx) , (6) 

where tx is the X-ray optical depth and F^. is: 

with L* being the luminosity of the central source. 

We estimate rx between the central source and a point in a flow from 

pr 

rx= Kxp(r, 9) dr, (8) 







where k,x is the absorption coefficient, and r is the distance from the central source. We 
assume kx = 0.4 g~^ cm^ which is numerically of value of electron scattering. We found that 
in most cases, the wind is optically thin and the wind column density. 



Nnie) = P^dT (9) 



is less than 10^^ cm ^. 



3. RESULTS 



To complete the specification of the simulations, we need to assign the properties of 
the disk atmosphere (the lower boundary condition for the wind solution). For the density. 
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we use a simple distribution at ^ = 90° of the form: p — Po('"/nc)~°, where a and po 
are constants. We run the model for a = 0, 1, 2, and 3 and for po = 10~^^, 10~^^, and 
10~^° g cm^^ (see table 1). The temperature is set to 10*^ K, while the radial velocity is set 
to zero. In addition, we enforce Keplerian rotation at = 90°. The size of the computational 
domain is defined in the following way: ri = 0.05 -Ric while Tq = 20i?ic. We follow M06 and 
NOG and adopt the following properties of the system: the mass of the central black hole Mbh 
= 7 Mq, the luminosity = 0.03 in units of the Eddington limit, LEdd; and the Compton 
temperature, Tic = 1-4 x 10^ K. The adopted corresponds to the mass accretion rate 
Ma = 4.4 X 10^^ g s-^ 

3.1. PROPERTIES OF THE FIDUCIAL RUN 

Our fiducial run (C8) has the following parameters: a=2, po = 10^^^ g cm^^, d6'k+i/(i^k ~ 
0.97, dr\j^i/dr\ = 1.04. This run settles to a steady state after about 3 sound-crossing times. 

We show the geometry and structure of the fiducial run in Figure 1. The top left 
panel shows a density distribution with a significant departure from spherical symmetry. 
The temperature (top right panel) is close to 0.1 Tic in the outer part of the wind. In a 
very narrow region around the rotational axis, the temperature is much less than Tic 
r > 16i?ic, whereas at the small radii the temperature is comparable to Tic- The X-ray 
irradiated gas accelerates rapidly and becomes supersonic relatively close to the disk. As 
shown in the bottom left panel, the contour for the Mach number of one corresponding to 
the wind launched at relatively large radii is almost a straight line making an angle of about 
15° with the disk midplane. Streamlines (bottom right panel) show that most of the disk 
wind follows almost perfect straight lines inclined at the angle of about 45° with respect to 
the disk midplane. The streamlines are curved only in zones near the rotational axis and 
very close to the disk. It appears that the outer wind is self-similar (we will return to this 
point at the end of this subsection). The streamlines also show that the flow expands in 
the angular direction, reminiscent of a spherical outflow, only for streamlines originating at 
r ^ 3 it!ic. 

An accretion disk is geometrically extended and as such it intercepts the central radiation 
over a large range of radii. However, the thermal expansion drives a signiflcant wind only 
within a relatively narrow radial range. One can show this wind property by plotting the 

product of the density and the velocity normal to the disk as a function of the radius along 
the disk midplane (see Figure 2). Most of the outflow in the fiducial run comes from the 
narrow ring on the disk in the zone between 0.07 i?ic and 0.2i?ic having its maximum at 
~ 0.1 Ric- For radii > 0.2 Ric the mass flux density scales like r"'', with q = 1.76. 
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The wind is launched in a non-spherical way and it remains non-spherical as it expands. 
To illustrate this point, Figure 3 shows various flow properties at as functions of 6. In 
particular, the top panel shows the product of density and radial velocity (solid line) and 
the accumulated mass- loss rate (dotted line). We evaluate the accumulative mass- loss rate 
throughout To using the following formula: 

re 

dm{e)=Anr^^ I pv^suiOdO. (10) 

One finds that most of the mass fiows out of the computational domain through the outer 
spherical boundary over a broad range of angles between 20° and 75°. The outflow for large 
9 contributes insignificantly. The total mass-loss rate for the fiducial run is M„ = (im(90°) = 
3.3 X lO^^gs^^. Note that this is larger than the assumed accretion rate by a factor of 
7.5. We will return to this point in Section 5. 

The middle panel of Figure 3 shows that the photoionization parameter (solid line) 
decreases with increasing 6^ which is another indication that the wind is not spherically 
symmetric. In particular, ^ changes from 10^ to 10^, whereas the column density (dotted 
line) increases from below 10^^ cm~^ to lO^^cm"^. 

The bottom panel of Figure 3 presents the radial velocity (solid line) and the number 
density (dotted line) vs. 9. It is evident from this panel why the zones close to the disk and 
the rotational axis contribute insignificantly to the total mass-loss rate: the zone near the 
disk is very dense but also very slow, while the zone near the axis of rotation has extremely 
low density, so that the product of the density and radial velocity is very small. 

To specify the departure of the disk wind from a purely radial wind. Figure 4 shows 
^ (top panel) and n (bottom panel) as a function of radius for various 9. In the case of 
an optically thin radial wind with a constant velocity, ^ is constant and n oc r^^. In our 
solution, ^ decreases with increasing radius especially for small and intermediate 9, whereas 
n decreases slower than r"^ and can even increase with radius for 9 < 62°. 

How does our solution compare to a /3- velocity law, in terms of approaching a constant 
value at large r? As expected, it is in agreement for large radii. Figure 5 shows an example 
of the radial velocity as a function of radius for 9 — 76" (this corresponds to the dot-dashed 
lines in Fig. 4). The radial velocity increases up to the Compton radius and then remains 
roughly constant. Also, there is a visible dip at 0.4 Ric, which is a consequence of the radial 
line crossing different streamlines. The function has a nearly flat shape for radii larger than 
l.Oi^ic. 

We conclude that a thermal disk wind can not be approximated by a radial outflow. Font 
et al. (2004) found that a disk wind expands quasi-spherically if a relatively steep decline 
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of the mass flux density with radius is assumed {pvo oc r"~^, with q <; 2.5). As we mentioned 
above, in our self-consistent simulations, q ~ 1.76 in the self-similar part of the wind. It is 
unfortunate that the wind expansion is non-radial because a radial outflow is a very simple 
case that can be (and has been) straightforwardly implemented into photoionization/spectra 
calculations. Is it then possible to approximate the thermal wind with some other simple 
model? It is beyond the scope of this paper to answer this question. 

3.2. COMPARISON OF THE FIDUCIAL RUN WITH OTHER RUNS 

To check the robustness of the solution presented above, we explored effects of the nu- 
merical resolution and the density along the lower boundary. We present the most significant 
properties of those runs in Table 1. The columns arc organized in the following way: (1) 
the name of the run; (2) the a parameter used in the density profile; (3) po, the normalized 
density at the lower boundary of the computational domain; (4) d9i^i/d9i, the angular reso- 
lution; (5) p(rin) the disk density at the inner radius; (6) the total wind mass-loss rate; 
(7) the maximal angle O^ax for which the integration of the total wind mass-loss is computed 
[Mw = dm{9max)]; (8) the maximal radial velocity on the outer shell; (9) the total time of 
the run. 

As mentioned in the Section 3, the density along the lower boundary {9 — 90°), is 
specified by two parameters, a and Po- We have run several models for various a and Po 
to check if our solution depends on the density along the disk. As for radiation driven disk 
winds (see e.g., Proga et al. 1998), we find that as long as the density along the disk is high 
enough, the gross properties of the thermal disk wind do not depend on the assumed disk 
density. 

For a relatively small density along the disk, the X-rays heat the gas to the Compton 
temperature even at = 90°. This means that the computational domain does not capture 
a cold disk and its part that is in hydrostatic equilibrium. Consequently, the simulations do 
not follow the transition between a cold and a hot fiow. Therefore, the mass flux density 
from the lower boundary is set not by the physics of a cold disk being heated by X-rays 
but by a choice of the density. Comparing M„ for runs C4, C8, Cll, and D2, we find that 
the mass-loss rate depends on the density if the density is too low. Only when the density 
along the disk is high enough that the central X-rays do not heat the gas along the lower 
boundary, the simulations capture the transition from the cold to hot phase. The mass-loss 
rate then becomes independent of a particular choice of the density in the cold part of the 
disk. 
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However, we also find that for very high densities along the disk there is a technical 
problem with a proper measurement of M^. Our goal is to compute a wind from a geo- 
metrically thin disk. But if we choose too a high density along the lower boundary, the 
disk will flare and its thickness at large radii becomes substantial. During the simulations, 
the dense disk remains cold and nearly in hydrostatic equilibrium but it will subsonically 
fluctuate. Therefore, when computing the wind mass-loss rate, one needs to keep in mind 
these fluctuations of the dense disk and exclude the region very close to the disk midplane. 
Otherwise the subsonic oscillations of the very dense disk material would be counted as the 
disk outflow and this will yield an erroneous estimate of the wind mass loss rate. In practice, 
when calculating the total outflow rate for some runs, we stop the integration in eq. 10 at 
^max < 90°, where the transition from the wind to the disk occurs. This transition is not 
difficult to identify as it occurs where the density sharply increases with 9 (see the bottom 
panel in Fig. 3). The seventh column in Table 1 lists the maximum angles we used. 

For the fiducial run as well as some other runs, the region near the disk midplane does 
not affect estimations of M„. This is one of the reasons we chose the run C8 with po — 
10~^^gcm^^ as the most suitable run in terms of having not too high a density, as only a 
small portion of the dense disk will enter the domain. On the other hand, the density cannot 
be so small as to make the transition from the disk to the wind outside the computational 
domain. All runs with the density comparable to or larger than that in the fiducial run have 
a total mass-loss rate of about 3.3 x lO^^gs"^. Other wind properties are also similar (see 
Table 1). 

As for the resolution study, we find that changes in the d6\^j^i / dOy, factor produced small 
differences of the final steady state and such effects just shghtly changed calculated physical 
quantities and properties of the solution. 

4. COMPARISON WITH OBSERVATIONS 

We start our comparison with the observations by considering the wind properties in- 
ferred from the observation and then we present an example of the synthetic spectrum 
directly compared with the observed X-ray spectrum. As we described in Section 1, the 
number density of the wind is inferred to be quite high, i.e., 5 x 10^^ cm~^. 

Figure 6 shows the scatter plot of the photoionization parameter vs. Vi- based on run 
C8. Dots correspond to n < 10^^ cm~^, whereas diamonds correspond n > 10^^ cm~^. Note 
that there are no diamonds corresponding to > lOOkms"^. This means that our model 
fails to predict the wind density and velocity as inferred from the observations. We stress 
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that in the computational domain of our simulation there is gas with densities as high as 
n ~ 5 X 10^^ cm~^, but this high density gas is very slow (e.g., the bottom panel of Fig. 3). 
This simple comparison shows that our results support M06 and Kallman et al.'s (2009) 
conclusion: a thermal wind can not account for the X-ray observations. 

As another test of our model, we computed synthetic spectra using the photoionization 
code XSTAR (Kallman &: Bautista 2001; Bautista & Kallman 2001). Our method of making 
this calculation is very similar to that used and described in Dorodnitsyn et al. (2008). We 
have computed spectra of the fiducial model for several inclination angles. We have also 
calculated the emissivity and opacity at each point in the hydrodynamic fiow at the end of 
the simulations, and then integrated the formal solution of the transfer equation along radial 
rays to get the spectrum. In doing this we ignore emission. Since emission will tend to fill in 
absorption lines, thus reducing their strength, these models represent an overestimate to the 
amount of absorption coming from such a fiow. For most inclinations, the spectra show no 
absorption features or only very weak ones. Figure 7 compares the count spectrum observed 
by the HETG during an X-ray bright phase of the GRO J1655-40 2005 outburst (see M06 
for more details) with an example of our synthetic spectrum. 

For this figure, we have chosen an almost 90° inchnation because it yields the best fit 
to the data and the right ionization for the most lines. Nonetheless, even for this very high 
inclination, the range of ionization is too narrow, and the ionization is probably lower than 
that observed. It is possible to obtain higher ionization for lower inclinations, but then the 
column density is lower, making the fit is even worse. To see if we can improve the fit for a 
given density structure, we have to put in the 600 km s~^ blueshift by hand while computing 
this spectrum (note that our wind model does not predict such a high velocity of the dense 
wind at this inclination angle). However, neither this nor any other simple model fits the Fe 
Ka fine profile showing the high velocity (~1600 km s~^) component. We found also that 
the curve-of- growth of the lines does not fit: there are no saturated lines in the model. Still, 
close to half of the observed fines are predicted in the model, including the Fe XXII fines. 

In addition to an inadequate mass loss rate, the model shows a geometric structure 
different from that required by spectra of GRO J1655-40. From the lack of emission lines, 
M06 and NOG concluded that the global covering factor of the wind is small, while Figure 3 
shows that the density is reasonably constant over a broad range of angles. 
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5. SUMMARY AND THE DISCUSSION 

In this paper, wc presented axisymmetric, time-dependent hydrodynamical numerical 
simulations of thermally driven disk winds. To apply our results to GRO J 1655-40, we 
adopted the observed mass of its black hole, and the properties of its X-ray radiation. We 
performed the simulations using the same code and in a similar fashion to PK02 who com- 
puted wind for low-mass X-ray binaries. The main difference is that we turned off the 
radiation driving here because as found by PK02 it is neghgible for the stellar black hole ac- 
cretors. Turning off the radiation driving makes the simulations run much faster. The main 
goal of our simulations is to check whether, in the absence of magnetic fields and within 
the given simplifications and assumptions, the thermally driven wind can account for the 
observation. 

To compare the results with the observations, we also computed transmitted X-ray 
spectra using the photoionization code XSTAR and our wind solution. Our main findings 

are: 

1) The density of the fast-moving part of the wind is more than one order of magnitude 
lower than that inferred from the observations. Consequently, contrary to the claim made 
by NOG, the thermal model fails to predict synthetic spectra with line absorption as strong 
and as blueshifted as those observed. Overall our results support the conclusion reached by 
Miller et al. (2008) and Kallman et al. (2009), that GRO J1655-40, and hkely other X-ray 
transient sources, have thermal winds insufficient to explain the observed spectra; 

2) Despite the thermal wind being weak and Gompton thin, the ratio between the mass-loss 
rate and the mass accretion rate is about seven; 

3) We discovered that the outer wind is self-similar. 

One should ask if our simulations are conclusive. Namely, can we rule out thermal 
driving as a mechanism responsible for the observed wind even though we did not consider 
some effects such as thermal conduction? Gonduction could increase the degree of disk 
heating and could in principle increase the gas density above the disk photosphere. However, 
it is unclear if this increase would lead to an increase in the density at the wind base and 
consequently in an increase of the density of the fast wind, which is what is needed to account 
for the observations. The density of the wind base also could be sensitive to the details of 
radiative transfer and spatial resolution. Hydrostatic models of the ionization layer on the 
surface of the disk by Jimenez-Garate et al. (2002, 2005) showed that there is much more 
opacity than the electron scattering opacity in the very thin layer that produces UV and 
most of the X-ray emission on top of the optically thick disk. Therefore, our pure electron 
scattering assumption might underestimate the opacity. However, if this were true then that 
would mean that thermal driving is even less efficient because the irradiation would penetrate 
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the disk even less. We expect that the mass loss rate is set at the sonic point, which would 
occur in the million Kelvin temperature plateau seen in the hydrostatic models. In that 
region, the models indicate that photoabsorption opacity due to Fe XXV and Fe XXVI is 
comparable to the electron scattering opacity. The radiative heating also can be affected by 
disk flaring or the shape of the disk, in general. For example, at a given radius a strongly 
flared disk might intercept more central radiation than a flat disk. Consequently, a resulting 
outflow might be more coUimated and denser. We intend to extend our model to explore 
these effects in near future. However, before that we plan to check if the inclusion of magnetic 
driving will produce a dense fast wind. 

It is important to remember that the outburst spectra of GRO J 1655-40 and GRS 
1915+05 are exceptional. Spectra of these objects in their low states, as well as spectra 
of other black hole binaries, show absorption lines of only Fe XXV and Fe XXVI. That 
suggests lower densities and mass loss rates, and it raises the question of whether thermally 
driven winds might account for those observations. The microquasar H1743— 322 (Miller et 
al. 2006) is a good example, with Fe XXV and Fe XXVI equivalent widths of up to 4.3 mA 
and 6.8 mA, respectively. The densities are below 10^^ cm^'^ and the ionization parameters 
are above 10^ in 4 separate observations. We have examined our models at sightlines farther 
from the disk, where the velocities match those reported by Miller et al. (2006). The Fe 
XXV and Fe XXVI column densities are too small to account for the observed equivalent 
widths by a factor of a few, and we conclude that thermally driven winds cannot account 
for the typical low state winds. 

In the models presented here, we found that M^/Ma ~ 7. This is a fairly large ratio, 
which could mean that the thermally driven wind can signiflcantly change the mass flow in 
the disk. Therefore, we decided to check if our result is consistent with the results of others 
who modeled thermally driven winds. Of particular interest is work by Woods et al. (1996), 
who studied thermally driven disk winds in great detail and performed many axisymmetric 
wind simulations for various luminosities. They summarized their results in an analytic 
formula which fitted the mass flux density distribution obtained from simulations for various 
(their eq. 5.2). 

To find M„ / Ma predicted by Woods et al.'s simulations, we integrated their formula 
over the disk surface to obtain Mwind and then expressed through the following equation 
of the accretion luminosity: 

GMe^ (11) 

where rg is the Schwarzschild radius, which yields a relation between M^j and M^. We found 
that Woods et al. simulations predict M„/ Ma in a range from 2.0 to 6.0 depending on the 
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luminosity. The L* dependence is weak for L*/LEdd < 0.1. Thus / Ma is higher than, 
but comparable, to Woods et al.'s results for L^, on the lower side of their luminosity range. 
Nevertheless, both set of simulations predict / > 1. 

Our conclusion from this analysis is that is greater than the rate at which the 
central engine is fueled. This in turn can make the disk variable. One could even expect 

this relatively strong wind to cause a recurrent disappearance of the inner disk even if the 
disk is fed at a constant rate at large radii. Neilsen & Lcc (2009) suggested such a process 
while interpreting jet /wind/radiation variability in the microquasar GRS 1915+105. It is 
beyond the scope of this work to model the effect of the wind on the disk. However, this 
problem was studied by Shields ct al. (1986), who showed that due to viscous processes no 
oscillations appear for / < 15! We conclude that the thermal wind is to weak to 
cause the oscillation. But the thermal wind is also too weak to account for the observed wind. 
Therefore, it is possible that in GRO J1655-40, and other sources, e.g., GRS 1915+105, the 
wind responsible for the observed X absorption will be so strong that / Ma > 15 and as 
such cause disk oscillations and contribute to the observed disk variability. 

We thank Tim Waters for his comments on the manuscript. We acknowledge support 
provided by the Chandra awards TM8-9004X and TMO-llOlOX issued by the Chandra X- 
Ray Observatory Center, which is operated by the Smithsonian Astrophysical Observatory 
for and on behalf of NASA under contract NAS 8-39073. 
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Table 1: Summary of the results: 



No. of run 


a 


Po 


d9i+i/d9i 


P {Tin} 






Vr 


*/ 






(gcm-3) 




(gcm-3) 


(iQiSgs-i) 


n 


(kms-1) 


(2.2 X lO^s *) 


A 





10-10 


0.97 


10-9 


0.33 


59 


587 


10 


B 


1 


10-10 


0.97 


10-9 


0.33 


67 


599 


10 


CI 


2 


10-12 


1.00 


10-10 


6.6 


90 


746 


10 


C2 


2 


10-12 


0.99 


10-10 


6.8 


90 


708 


10 


C3 


2 


10-12 


0.98 


10-10 


6.9 


90 


675 


10 


C4 


2 


10-12 


0.97 


10-10 


7.0 


90 


642 


10 


C5 


2 


10-12 


0.96 


10-10 


7.1 


90 


619 


10 


C6 


2 


10-12 


0.95 


10-10 


7.1 


90 


604 


30 


C7 


2 


10-11 


0.99 


10-9 


0.32 


90 


703 


10 


C8 


2 


10-11 


0.97 


10-9 


0.33 


90 


627 


10 


C9 


2 


10-11 


0.95 


10-9 


0.33 


90 


587 


10 


CIO 


2 


10-10 


0.99 


10-8 


0.33 


63 


692 


30 


Cll 


2 


10-10 


0.97 


io-« 


0.33 


62 


603 


30 


C12 


2 


10-10 


0.95 


10-8 


0.33 


63 


553 


10 


Dl 


3 


10-11 


0.97 


10-8 


0.11 


90 


670 


10 


D2 


3 


10-10 


0.97 


10-7 


0.32 


90 


621 


30 



* 2.2 X 10® s is the sound crossing time = (ro/cs), calculated for T=ric. 
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Fig. 1. — The fiducial run C8. Top left panel: Color density map. Top right panel: 
Color temperature map. Bottom left panel: The Mach number M contours (the poloidal 
component only). Contours are for M = 1, 2, and 3 (bottom to top). Bottom right panel: 
The flow streamlines. Note the self-similarity of the streamlines especially those arising from 
the outer disk. In all panels the rotation axis of the disk is along the left-hand vertical frame, 
while the midplane of the disk is along the lower horizontal frame. Lengths are expressed in 
units of the Compton radius, Ric- 
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Fig. 2. — The mass flux density, pve^ as the function of radius along the disk midplane, 
e = 90° for the fiducial run. 
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Fig. 3. — Quantities at the outer boundary, To of the fiducial model. The ordinate on the 
left hand side refers to the sohd fine, while the ordinate on the right hand side refers to 
the dotted line. Top panel: Mass flux density and the accumulative mass-loss rate. Middle 
panel: Photoionization parameter and the column density. Bottom panel: radial velocity 
and number density. 
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Fig. 4. — The radial profiles of the photoionization parameter (top panel) and the number 
density (bottom panel) for seven polar angles along the radius for the fiducial run. The angles 
are: 9 = 48.3° (thick soUd), 6 = 60.5° (dotted), 6 = 69.4° (dashed), 6 = 76.0° (dot-dashed), 
e = 80.9° (triple dot-dashed), 6 = 84.5° (long dashed) and 9 = 89.1° (thin sohd). Note that 
a simple spherically expanding outflow poorly approximates the disk wind solution (see the 
main text for more details). 
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Fig. 5. — Radial velocity as the function of radius at ^ = 76° for the flducial run. The 
non-monotonic radial profile demonstrates that the disk wind is not radial (see also Fig. 4). 
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Fig. 6. — Scatter plot of photoionization parameter vs. v.,. for the fiducial run. Dots 
correspond to n < 10^^ cm~^ whereas diamonds correspond n > 10^^ cm~^. The lack of 
points corresponding to Vr > 100 km s~l means that the thermal wind can not account for 
a high density fast outflow as observed in GRO J1655-40 during an X-ray-bright state. 
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